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We present a measurement of direct photon pair production cross sections using 4.2 fb _1 of data 
collected with the DO detector at the Fermilab Tevatron pp Collider. We measure single differential 
cross sections as a function of the diphoton mass, the transverse momentum of the diphoton system, 
the azimuthal angle between the photons, and the polar scattering angle of the photons. In addition, 
we measure double differential cross sections considering the last three kinematic variables in three 
diphoton mass bins. The results are compared with different perturbative QCD predictions and 
event generators. 

PACS numbers: 13.85.Qk, 12.38.Qk 



At a hadron collider, the direct photon pair (DPP) pro- 
duction with large diphoton invariant mass (M 77 ) cons- 
titutes a large and irreducible background to searches 
for the Higgs boson decaying into a pair of photons, 
for both the Fermilab Tevatron Q and the CERN LHC 
experiments 0. DPP production is also a significant 
background in searches for new phenomena, such as new 
heavy resonances Q , extra spatial dimensions [ij , or cas- 
cade decays of heavy new particles Thus, precise mea- 
surements of the diphoton differential production cross 



sections for various kinematic variables and their theo- 
retical understanding are extremely important for future 
Higgs and new phenomena searches. 

In addition, DPP production is interesting in its own 
right, and is used to check the validity of the predictions 
of perturbative quantum chromodynamics (pQCD) and 
soft-gluon resummation methods implemented in theo- 
retical calculations. Measurements involving the dipho- 
ton final state have been previously carried out at fixed- 
target [fj, |7[ and collider 8 10] experiments. However, 
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the large integrated luminosity accumulated by the DO 
experiment in pp collisions at yfs = 1.96 TeV at the 
Fermilab Tevatron Collider allows us to perform precise 
measurements of several observables in kinematic regions 
previously unexplored, as well as, for the first time, the 
measurement of double differential cross sections for this 
process. 

The DPP events produced in pp — > 77+X are expected 
to be dominantly produced via qq scattering (qq — !• 77) 
and gluon-gluon fusion (gg — > 77) through a quark-loop 
diagram. In spite of the suppression factor of a 2 s for 
99 ~~ ^ 77 as compared to qq — > 77, the former still gives a 
significant contribution in kinematic regions where the gg 
parton luminosity is high, especially at low M 77 . Figure[T] 
shows the expected contribution to the total DPP rate 
from gg — > 77, as predicted by the PYTHIA ll| Monte 
Carlo (MC) event generator with the CTEQ6.1L parton 
distribution function (PDF) set 12]. In addition, direct 
photons may result from single or double fragmentation 
processes of the partons produced in the hard scattering 
[T3L Hij ]. However, a strict photon isolation requirement 
significantly reduces the rate for these processes. 




50 100 150 200 250 300 350 

M„ (GeV) 

FIG. 1: The fraction of events produced via gg — > 77 scat- 
tering relative to total diphoton production as a function of 
M 77 , as predicted by the pythia event generator using the 
CTEQ6.1L PDF set. Photons are required to have transverse 
momentum pr > 21(20) GeV for the highest (next-to-highest) 
Pt photon and pseudorapidity \rj\ < 0.9 [l5| |. 



In this Letter, we present measurements of the DPP 
production cross sections using data collected by the DO 
experiment from August 2006 to June 2009. The cross 
sections are measured differentially as a function of M 77 , 
the diphoton transverse momentum (pj 7 ), the azimuthal 
angle between the photons (A(/> 77 ), and the cosine of 
the polar scattering angle of the photon in the frame 
with no net transverse momentum of the diphoton sys- 
tem (defined as cos 6** = tanh^i — 77 2 )/2] , where r/i^) 
is the pseudorapidity of the highest (next-to-highest) pr 
photon). These kinematic variables probe different as- 
pects of the DPP production mechanism. For instance, 
the shapes of the pl^ and A0 77 distributions are mostly 



affected by initial state gluon radiation and fragmenta- 
tion effects. In addition, the M 77 spectrum is particu- 
larly sensitive to potential contributions from new phe- 
nomena. The cos 9* distribution probes PDF effects and 
the angular momentum of the final state, which should 
be different for QCD-mediated production as compared, 
for example, to the decay of a spin-0 Higgs boson fDij . 
The measured cross sections__are compared to theoretical 
predictions from RESBOS 



14(, and pythia 



[13[, DIPHOX 

[ill ]. Both RESBOS and diphox provide next-to- leading 
order (NLO) predictions in pQCD, however the gg — > 77 
contribution is considered only at leading order (LO) in 
diphox. PYTHIA is a parton shower MC event generator 
that includes the above processes at LO. In diphox, the 
explicit parton-to-photon fragmentation functions are in- 
cluded at NLO, while in RESBOS a function approximat- 
ing rate from the NLO fragmentation diagrams is intro- 
duced. Also, only in RESBOS, the effects of soft and 
collinear initial state gluon emissions are resummed to 
all orders. This is particularly important for the descrip- 
tion of the pip (A0 77 ) distribution, which is a (S-function 
at LO and diverges at NLO as pl^ — > (A0 77 — > w). 

The DO detector is a general purpose detector dis- 
cussed in detail elsewhere [lj|. The subdetectors most 
relevant to this analysis are the central tracking system, 
composed of a silicon microstrip tracker (SMT) and a 
central fiber tracker (CFT) embedded in a 2 T solenoidal 
magnetic field, the central preshower detector (CPS), and 
the calorimeter. The CPS is located immediately before 
the inner layer of the calorimeter and is formed of approx- 
imately one radiation length of lead absorber followed by 
three layers of scintillating strips. The calorimeter con- 
sists of a central section with coverage in pseudorapidity 
of |??dot| < 1.1 15|, and two end calorimeters covering up 



to |r?det| ~ 4.2. The electromagnetic (EM) section of the 
calorimeter is segmented longitudinally into four layers 
(EMz, i — 1,4), with transverse segmentation into cells 
of size A?7dct x A^det = 0.1 x 0.1 [15| , except EM3 (near 
the EM shower maximum), where it is 0.05 x 0.05. The 
calorimeter is well-suited for a precise measurement of 
the energy and direction of electrons and photons, pro- 
viding an energy resolution of about 3.6% at an energy 
of 50 GeV and an angular resolution of about 0.01 radi- 
ans. The energy response of the calorimeter to photons 
is calibrated using electrons from Z boson decays. Since 
electrons and photons shower differently in matter, addi- 
tional corrections as a function of r\ are derived using a 
detailed GEANT-based 17f simulation of the DO detector 
response. These corrections are largest [(2.0 — 2.5)%] at 
low photon energies (« 20 GeV). The data used in this 
analysis were collected using a combination of triggers re- 
quiring at least two clusters of energy in the EM calorime- 
ter with loose shower shape requirements and varying pr 
thresholds between 15 GeV and 25 GeV, and correspond 
to an integrated luminosity of 4.2 ± 0.3 fb _1 [l8j |. 

Events are selected by requiring two photon candi- 
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dates with transverse momentum pr > 21 (20) GeV for 
the highest (next-to-highest) pt photon candidate and 
pseudorapidity \rj\ < 0.9, for which the trigger require- 
ments are > 96% efficient. The minimum p? require- 
ments for the two photon candidates are chosen to be 
different following theoretical discussions IH 14 1 and a 
previous measurement [l(| • The photon px is computed 
with respect to the reconstructed event primary vertex 
(PV) with the highest number of associated tracks. The 
PV is required to be within 60 cm of the center of the de- 
tector along the beam axis. The PV has a reconstruction 
efficiency of about 98% and has about 65% probability 
of being the correct vertex corresponding to the hard 
PP 77 + X production. 

Photon candidates are formed from clusters of 
calorimeter cells within a cone of radius 1Z = 
^{Arj) 2 + (A0) 2 = 0.4 around a seed tower [16]. The 
final cluster energy is then recalculated from the inner 
core with 7Z — 0.2. The photon candidates are selected 
by requiring: (i) > 97% of the cluster energy be deposited 
in the EM calorimeter layers; (ii) the calorimeter isola- 
tion X = [E tot (0A) -£ em (0.2)]/£em(0.2) < 0.10, where 
Etat(7Z) [Eem(T^)] is the total [EM only] energy in a cone 
of radius 1Z; (iii) the px scalar sum of all tracks origi- 
nating from the PV in an annulus of 0.05 < 1Z < 0.4 
around the EM cluster be < 1.5 GeV; and (iv) the energy- 
weighted EM shower width be consistent with that ex- 
pected for an electromagnetic shower. To suppress elec- 
trons misidentified as photons, the EM clusters are re- 
quired to not be spatially matched to significant tracker 
activity, either a reconstructed track or a density of hits 
in the SMT and CFT consistent with that of an elec- 
tron [19(. In the following, this requirement will be re- 
ferred to as the "track-match veto" . 

To further suppress jets misidentified as photons, an 
artificial neural network (NN) discriminant which ex- 
ploits differences in tracker activity and energy deposits 
in the calorimeter and in the CPS between photons and 
jets is defined [l|. The NN is trained using 7 and jet 
PYTHIA MC samples. The shapes of the NN output 
(Onn), normalized to unit area and obtained after ap- 
plying all data selection criteria, are shown in Figure 
[2j exhibiting a significant discrimination between pho- 
tons and jets. Photon candidates satisfy the requirement 
Onn > 0.3, which is «98% efficient for photons and re- 
jects s» 40% of the jets misidentified as photons. The 
Onn shape is validated in data. For photons a data sam- 
ple consisting of photons radiated from charged leptons 
in Z boson decays (Z — » £ + £~j, £ = e,fj) [20j is used. 
The MC modeling of the Onn shape for jets is validated 
in a sample of photon candidates selected by inverting 
the photon isolation (T > 0.07), a requirement that sig- 
nificantly enriches the sample in jets. The data and MC 
Onn shapes are compared in Figures [2] and [3] and found 
to be in good agreement. 

Finally, the two photon candidates are required to be 
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FIG. 2: Comparison of the normalized Onn spectra for pho- 
tons from DPP MC and Z —¥ £ + £~y data and for misidentified 
jets from dijet MC. 



spatially separated from each other by a distance in r\ — <fi 
space AT? > 0.4 and to satisfy M 77 > pip . The lat- 
ter requirement is satisfied by the majority (« 92%) of 
DPP events and, together with the photon isolation re- 
quirements, allows significant suppression of the contri- 
bution from the fragmentation diagrams, thus restricting 
the data-to-theory comparison to the region where the 
theoretical calculations should have smaller uncertainties 

After imposing all requirements, 10938 events with 
diphoton candidates are selected in data. This sam- 
ple includes instrumental background contributions from 
7+jet and dijet production, where a jet is misidentified 
as a single photon as a result of fluctuations in the parton 
fragmentation into a well-isolated neutral meson (ir° or rf) 
decaying into a final state with two or more photons. An 
additional smaller background contribution results from 
Z-boson/Drell-Yan production events Z/-f* —> e + e~ 
(ZDY) in which both electrons are misidentified as pho- 
tons. 

The contribution from ZDY events is estimated us- 
ing the MC simulation with pythia, normalized to the 
NNLO cross section 2l|. The selection efficiencies de- 



termined from the MC simulation are corrected to those 
measured in the data. On average, each electron has 
a 2% probability of satisfying the photon selection cri- 
teria, mainly due to the inefficiency of the track-match 
veto requirements. The total ZDY contribution is es- 
timated to be 161 ± 20 events. Backgrounds due to 
7+jet and dijet events are estimated from data by using 
a 4 x 4 matrix background estimation method [l(. Af- 
ter applying all of the selection criteria described above, 
a tighter Onn requirement (Onn > 0.6) is used to 
classify the data events into four categories, depend- 
ing on whether both photon candidates, only the high- 
est pt one, only the next-to-highest px one, or nei- 
ther of the two photon candidates pass (p) or fail (/) 
this requirement. The corresponding number of events 
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(after subtraction of the estimated ZDY contribution) 
compose a 4-component vector (N pp ,Npf,Nfp,Nff). 
The difference in relative efficiencies of the NN > 
0.6 requirement between photons and jets allows es- 
timation of the sample composition by solving a lin- 
ear system of equations: (N pp ,N p f,Nf p ,Nff) T = £ x 
(-/V 77 , Njj, iVj 7 , Njj) T , where iV 77 (Njj) is the number 
of DPP (dijet) events and iV 7 j (Nj 7 ) is the number of 
7+jet events with the (next-to-)highest px photon can- 
didate being a photon. The 4x4 matrix £ contains 
the photon e 7 and jet £j e t efficiencies, estimated using 
photon and jet MC samples and validated in data. The 
efficiencies are parameterized as a function of the pho- 
ton candidate r\ and vary within (90 — 95)% for e 7 and 
within (66 — 70)% for £j et . The systematic uncertainty 
on e 7 is estimated to be f .5% from a comparison of the 
efficiency as a function of rj between data and MC using 
samples of electrons from Z boson decays and photons 
from radiative Z boson decays. In order to estimate the 
systematic uncertainty on £j et , two independent control 
data samples enriched in jets misidentified as photons are 
selected, either by inverting the photon isolation variable 
(I > 0.07), or by requiring at least one track in a cone 
of 7Z < 0.05 around the photon, while keeping the re- 
maining photon selection criteria unchanged. In both 
cases the agreement with the MC prediction for £j et is 
found to be within 10%, which is taken as the system- 
atic uncertainty. The total number of DPP events is 
found to be iV 77 = 7307 ± 3I2(stat.), corresponding to 
an average DPP purity of ss 67%. Following this pro- 
cedure, the number of DPP events is estimated in each 
bin of the four kinematic variables considered (M 77 , pip , 
A0 77 , and |cos#*|). The largest kinematic dependence 
of the DPP purity is in terms of M 77 , with a variation 



between « 60% at 



40 GeV and close to 100% for 



M 77 > 200 GeV. As a function of the other kinematic 
variables, the DPP purity varies in the (60 — 70)% range. 
The relative systematic uncertainty on the purity results 
from the systematic uncertainties on e 7 and £j e t, and typ- 
ically varies within (II — 15)%. As a cross-check, the DPP 
purity was also estimated via a fit to the two-dimensional 
distribution in data of Onn, 71 versus Onn, 72 using tem- 
plates constructed from photons and jets in MC. The 
result was found to be in good agreement with that from 
the 4x4 matrix method. 

The estimated number of DPP events per bin is cor- 
rected for the DPP event selection efficiency and accep- 
tance. The selection efficiency is calculated using DPP 
events generated with PYTHIA and processed through a 
GEANT-based simulation of the DO detector. In order to 
accurately model the effects of multiple pp interactions 
and detector noise, data events from random pp cross- 
ings with a similar instantaneous luminosity spectrum as 
considered in the data analysis are overlaid on the MC 
events. These MC events are then processed using the 
same reconstruction code as for the data. Small differ- 
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FIG. 3: Comparison of the normalized Onn spectra for jets 
misidentified as photons in data and in dijet MC. 



ences between data and MC in the per-photon selection 
efficiencies are corrected for with suitable scale factors 
derived using control samples of electrons from Z boson 
decays, as well as photons from the radiative Z boson de- 
cays. The overall DPP selection efficiency after applying 
all selection criteria is estimated as a function of the vari- 
able of interest. In the case of pip ', A0 77 , and | cos#*|, 
it is about 64% with a (2 - 3)% variation across the bins, 
while for M 77 , the efficiency grows from about 60% at 
30 < M 77 < 50 GeV to 69% at M 77 > 200 GeV. The to- 
tal relative systematic uncertainty on the DPP selection 
efficiency is 4.3%, dominated by the track-match veto and 
photon Onn selections. The acceptance is calculated us- 
ing DPP events generated with RESBOS and is driven by 
the selections in rj^et (|?7det| < 0.9, applied to avoid edge 
effects in the central calorimeter region used for the mea- 
surement) and 4>det (to avoid periodic calorimeter module 
boundaries [16j | that bias the EM cluster energy and po- 
sition measurements), PV misidentification, photon en- 
ergy scale, and bin-to-bin migration effects due to the 
finite energy and angular resolution of the EM calorime- 
ter. The overall DPP acceptance varies within (45 — 64)% 
with a relative systematic uncertainty of (4 — 7)%. 

The differential cross sections da/dM 17 , dojdpip , 
da/dA4>^j, and da/d\ cos#*| are obtained from the num- 
ber of data events corrected for the background contribu- 
tion, divided by the trigger, vertex and diphoton selec- 
tion efficiencies, acceptance, integrated luminosity, and 
the bin width for each kinematic variable. The measured 
differential cross sections, compared to the theoretical 
predictions from RESBOS, are presented in Table Q] The 
average value for each variable in a bin was estimated us- 
ing RESBOS. The statistical uncertainty 5 sta t corresponds 
to the statistical precision on 7V 77 estimated in the 4x4 
matrix method, which can be sizable when values of e 7 
and ej e t are numerically close. 

Figure |4] shows a comparison of the measured differ- 
ential cross sections to the theoretical predictions from 
RESBOS, diphox, and PYTHIA. Systematic uncertainties 
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in the measured cross sections have large (> 90%) bin-to- 
bin correlations. There is a common 7.4% normalization 
uncertainty, resulting from the photon selection criteria 
(4.3%) and luminosity measurement (6.1%), that is not 
shown on the data points. The predictions from RES- 
BOS and DIPHOX are computed using the CTEQ6.6M 
PDF set 12], the DSS set of fragmentation functions 
22j, and setting renormalization jiR, factorization [ip, 
and fragmentation /j,f scales as /ir = \ip = A*/ = M 77 . 
The uncertainty due to the scale choice is estimated by 
simultaneous variation by a factor of two of all scales 
relative to the default choice and found to be about 
10% for M 77 and |cos6>*| and up to (15 - 20)% for 
high plf and low A</> 77 . The PDF uncertainty is esti- 
mated using DIPHOX and the 44 eigenvectors provided 
with the CTEQ6.6M PDF set [pf and found to be 
within (3 — 6)% for all four cross sections. The pre- 
dictions from PYTHIA are computed with "Tune A" llj . 
which uses the CTEQ5L PDF set. All theoretical predic- 
tions are obtained using diphoton event selection criteria 
equivalent to those applied in the experimental analy- 
sis. In particular, the photon isolation is required to be 
E™° = Eg* (OA) - £ 7 < 2.5 GeV, where £^ ot (0.4) is the 
total transverse energy within a cone of radius 1Z = 0.4 
centered on the photon, and Ej, is the photon transverse 
energy. For resbos and diphox, Etf* is computed at the 
parton level, whereas in the case of PYTHIA, it is com- 
puted at the particle level. This requirement suppresses 
the contributions from photons produced in the fragmen- 
tation processes and leads to a more consistent compari- 
son with the experimental result. Studies performed us- 
ing DIPHOX indicate that the contribution to the overall 
cross section from one- and two-fragmentation processes 
does not exceed 16% and significantly drops at large M 77 , 
p7p and small A0 77 to (1-3)%. In order to allow a direct 
comparison to the data, the NLO QCD cross sections 
obtained with RESBOS and DIPHOX are further corrected 
for contributions from multiple parton interactions and 
hadronization, both of which affect the efficiency of the 
isolation requirement. These corrections are estimated 
using DPP events simulated in pythia using Tunes A 
and SO The corrections vary within (4.0 — 5.5)% as 
a function of the measured kinematic variables and are 
consistent for both tunes within 0.5%. 

The results obtained show that none of the theoretical 
predictions considered is able to describe the data well 
in all kinematic regions of the four variables, resbos 
shows the best agreement with data, although systematic 
discrepancies are observed at low M 77 , high p 7 -, 7 , and 
low A(/> 77 . However, the agreement between RESBOS and 
data is fair at intermediate M 77 (50 — 80 GeV) , and good 
at high M 77 (> 80 GeV). The large discrepancy between 
RESBOS and diphox in some regions of the phase space is 
due to absence of all-order soft-gluon resummation and 
accounting gg — > 77 contribution just at LO in diphox. 

Further insight on the dependence of the pj. 7 , A<^ 77 , 



TABLE I: The measured differential cross sections in bins of 



M T) 77 

iw 77 , p T , 



A</> 77 , and I cos ( 



The columns 5 sta t and 8, 



syst 



represent the statistical and systematic uncertainties, respec- 
tively. Also shown are the predictions from RESBOS. 



Mv 



<M 7 



dcr/dM 77 (pb/GeV) 



(GeV) 


(GeV) 


Data 


#stat (%) (Ssyst 


(%) 
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30 - 
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3.11x10" 
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+26/ 


-29 
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50 


47.6 


1.74x10" 


1 
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-19 


1.22x10" 
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1.19x10" 
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2 


10 


+17/ 


-15 


4.09x10" 
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20 
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3 
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4 
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-14 
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4 
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da/dpf (pb/GeV) 



(GeV) 


(GeV) 


Data 


<5stat (%) <5syst 


(%) 


RESBOS 




0.0 - 


2.5 


1.5 


1.92x10" 


1 


15 


+18/ 


-19 


2.63x10" 
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5.0 
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-17 
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2.99x10" 
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13 
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-17 
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2 
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A<^> 77 
(rad) 



(Acp-f-y) d<j/dA<?!) 77 (pb/rad) 

(rad) Data <5 sta t (%) S sy st (%) resbos 
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- 1.88 


1.75 


4.32xl0 _1 
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I cos 6* (| cos #*|) 



da/d\ cos 9* |(pb) 
Data 5 Btat (%) (? syst (%) RESBOS 



0.0 


- 0.1 


0.05 


13.8 


8 
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17 


9.22 


0.1 
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9 


+17/- 


16 
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0.2 
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0.25 
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10 


+18/- 


16 
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0.3 


- 0.4 


0.35 


6.38 


12 


+17/- 


16 


5.90 


0.4 


- 0.5 


0.45 


4.77 


14 


+17/- 


16 


4.54 


0.5 


- 0.7 


0.57 


2.35 


15 


+17/- 


16 


2.16 



and I cos 6* | kinematic distributions on the mass scale 
can be gained through the measurement of double dif- 
ferential cross sections. For this purpose, the differential 
cross sections as functions of pj. 7 , A</> 77 , and | cos#*| are 

50 GeV, 50 - 80 GeV 
350 GeV. The results are presented in Tables 



measured in three M 77 bins: 30 



and 80 



corresponding to each of the three M 77 inter- 
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FIG. 4: The measured differential diphoton production cross sections as functions of (a) M 77 , (b) pip , (c) A</> 77 , and (d) | cos#*|. 
The data are compared to the theoretical predictions from RESBOS, diphox, and pythia. The predictions from RESBOS, and 
diphox use the CTEQ6.6M PDF set [12] and renormalization, factorization, and fragmentation scales /j,r — [if = \if = M 77 , 
while pythia uses the Tune A settings. Theoretical predictions are obtained using the following selections: two photons with 



p T > 21(20) GeV, |t?| < 0.9, 30 



<M 77 < 350 GeV, 



M 77 > pip , 



All > 0.4, A0 77 > 0.5tt, and < 2.5 GeV. The ratio of 



differential cross sections between data and RESBOS are displayed as black points with uncertainties in the bottom plots. The 
inner line for the uncertainties in data points shows the statistical uncertainty, while the outer line shows the total (statistical 
and systematic added in quadrature) uncertainty after removing the 7.4% normalization uncertainty. The solid (dashed) line 
shows the ratio of the predictions from diphox (pythia) to those from RESBOS. In the bottom plots, the scale uncertainties 
are shown by dash-dotted lines and the PDF uncertainties by shaded regions. 



vals. Each table is split into three sub-tables, showing 
results separately for d 2 a j dM 11 dp l p i , d 2 a/dM 11 dA(f) 11 , 
and d 2 a / dM 11 d\ cos 9* |. The measured cross sections 
for the pip, A</> 77 , and | cos 8*\ variables in the three 
mass bins are shown in Figures [5] — [7] and compared to 
the theoretical predictions. These results confirm that 
the largest discrepancies between data and RESBOS for 
each of the kinematic variables originate from the lowest 
M 11 region (M 77 < 50 GeV) . As shown in Figure [TJ this 
is the region where the contribution from gg — > 77 is 
expected to be largest. The discrepancies between data 
and RESBOS are reduced in the intermediate Af 77 region 
(50 — 80 GeV), and a quite satisfactory description of 



all kinematic variables is achieved for the M 77 > 80 GeV 
region, the relevant region for the Higgs boson and new 
phenomena searches. However, it should be pointed out 
that at the Tevatron, DPP production at high masses is 
strongly dominated by qq annihilation, in contrast with 
the LHC, where the contribution from gg and qg initiated 
process will be significant. It remains to be seen whether 
the addition of NNLO corrections to RESBOS, as done in 
(23} . will improve the description of the high pip (low 
A^> 77 ) spectrum at low M 77 . 

In summary, we have presented measurements of sin- 
gle and double differential cross sections for DPP pro- 
duction in pp collisions at y/s = 1.96 TeV. This analysis 



9 



D0, 4.2 fb" 1 

30 GeV ■ 




(a) 

M„ < 50 GeV 
data 

RESBOS 
DIPHOX 
PYTHIA 
PDF uncert. 
scale uncert. 



Sio- 1 



-i 

■D 



5 



10" 



O 3 

m 25 

uj 2 

°= 1.5 

S 1 

o 0.5 
« 



D0, 4.2 ft) -1 (b) 
30 GeV < < 50 GeV 



1.6 1.8 2 2.2 2.4 




o10 1 



■a 



10" 

2.5 
2 

1.5 
1 

0.5 



D0, 4.2 fb" 1 
30 GeV < M„ . 



(c) 



50 GeV 




} I 



FIG. 5: 



10 15 20 25 30 35 40 45 50 

P^tGeV) 

The measured double differential diphoton production cross sections as functions of 



0.1 0.2 0.3 0.4 0.5 



0.6 0.7 

I cose* | 



(a) 



pV 
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for 50 < M 77 < 80 GeV. The notations for points, lines and shaded regions are the same as in Figure 2] 
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uses 4.2 fb^ 1 of DO data, representing a twenty-fold in- 
crease in statistics relative to the last published Tevatron 
results jiol ]. The measured cross sections are compared 
to predictions from RESBOS, DIPHOX and PYTHIA, show- 
ing the necessity of including higher order corrections 
beyond NLO as well as the resummation to all orders of 
soft and collinear initial state gluons. These results allow 
the tuning of the theoretical predictions for this process, 
which is of great relevance for improving the sensitivity of 
searches for the Higgs boson and other new phenomena 
at the Tevatron and the LHC. 



TABLE III: The measured double differential cross sections 
in bins of , A</> 77 , and | cos#*|, in the region 50 < M 77 < 
80 GeV. The notations are the same as in Table ITU 
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TABLE II: The measured double differential cross sections in 
bins of pj 7 , A0 77 , and |cos6**|, in the region 30 < M 77 < 
50 GeV. The columns 5 s t a t and 5 aya t represent the statistical 



and systematic uncertainties, respectively, 
the predictions from RESBOS. 
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TABLE IV: The measured double differential cross sections 
in bins of p^f , A</> 77 , and | cos#*|, in the region 80 < M 77 < 
350 GeV. The notations are the same as in Table HT1 
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